Robust identification of isotropic diffuse gamma rays from Galactic dark matter 
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Dark matter annihilation in Galactic substructure will produce diffuse gamma-ray emission of 
remarkably constant intensity across the sky, making it difficult to disentangle this Galactic dark 
matter signal from the extragalactic gamma-ray background. We show that if Galactic dark matter 
contributes a modest fraction of the measured emission in an energy range accessible to the Fermi 
Gamma-ray Space Telescope, the energy dependence of the angular power spectrum of the total 
measured emission could be used to confidently identify gamma rays from Galactic dark matter 
substructure. 
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Introduction. — Many independent observations in- 
dicate the presence of a substantial non-baryonic dark 
matter component in the universe, but its nature remains 
unknown. Among the favored candidates for cold dark 
matter are weakly interacting massive particles such as 
the neutralino in supersymmetric theories and the light- 
est Kaluza-Klein particle in theories of universal extra 
dimensions. In these scenarios dark matter can annihi- 
late, producing standard model particles, and therefore 
may be detectable indirectly via its annihilation prod- 
ucts, such as gamma rays. The Large Area Telescope 
aboard the recently launched Fermi Gamma-ray Space 
Telescope represents an unprecedented leap in gamma- 
ray detection capabilities in the GcV energy band. Fermi 
will provide for the first time a detailed all-sky measure- 
ment of diffuse emission from ~ 20 GeV up to several 
hundred GeV, an energy range which is of particular in- 
terest for indirect dark matter searches. 

Gamma rays from Galactic dark matter are produced 
by annihilation in the smooth dark matter halo and in 
its substructure. While a few individual subhalos may be 
detectable as point sources, most will not be resolved by 
Fermi and will instead contribute to the measured diffuse 
emission [l|, |2j . Numerical simulations have shown that 
the radial distribution of subhalos is not concentrated 
toward the Galactic Center, but instead is quite extended 
relative to the smooth dark matter halo distribution 0, 0, 
[E00]- This leads to the important result that the diffuse 
emission from annihilation in Galactic substructure as 
observed from the Earth appears nearly isotropic on large 
angular scales @, d, Q . The diffuse flux from substructure 
typically exceeds that from the smooth halo for angles > 
a few tens of degrees from the Galactic Center [l[ 0, 0, 
[lo| . so substructure is the dominant contributor to the 
Galactic dark matter flux over most of the sky. 

Since gamma rays from Galactic substructure would 
appear in Fermi data as part of the isotropic diffuse 
emission, it is essential to identify ways of separating 
this signal from the extragalactic gamma-ray background 



(EGRB). Although the dark matter energy spectrum is 
distinct from that of known extragalactic source classes, 
confident identification of a dark matter signal on spec- 
tral grounds alone is unlikely. At energies as low as tens 
of GeV, photon-photon interactions with the UV, optical, 
and IR backgrounds (the extragalactic background light, 
EBL) lead to a suppression of the measured extragalactic 
source spectrum which can mimic the exponential cut-off 
in the dark matter spectrum which occurs at the energy 
corresponding to the dark matter particle mass m x . This 
problem is especially severe because the shape and char- 
acteristic energies of these spectral features are minimally 
constrained due to uncertainties in the EBL as well as 
in the properties of extragalactic source populations and 
dark matter. 

Recently, several studies have considered anisotropies 
in the diffuse gamma-ray background as a means of 
distinguishing between plausible source classes. The 
angular power spectrum of the diffuse emission has 
been predicted for various extragalactic source classes 
[lfjL 11, 12, 13 jl4. ljf a s well as for Galactic dark matter 



substructure 

In this Letter, we outline a general strategy for identi- 
fying contributions from multiple source populations to 
the diffuse gamma-ray emission using the energy depen- 
dence of the angular power spectrum, and consider here 
the specific challenge of extracting a Galactic dark mat- 
ter signal from the EGRB. Recognizing that (1) diffuse 
emission from unresolved extragalactic sources has differ- 
ent small-scale anisotropy properties than emission from 
Galactic dark matter substructure, and (2) their relative 
contributions to the total isotropic diffuse emission vary 
with energy, we show that the contribution from Galactic 
dark matter substructure to the total measured emission 
is identifiable as a modulation of the total angular power 
spectrum as a function of energy. 

The intensity energy spectrum. — The intensity 
energy spectrum Ie describes the photon intensity (pho- 
tons per area per time per solid angle per energy) as a 
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function of energy E. 

For dark matter annihilation, we consider the contin- 
uum gamma-ray intensity energy spectrum for which we 
use the analytic approximation for neutralinos given in 
[13] ■ The gamma-ray flux from annihilation is propor- 
tional to the dark matter density squared, so the nor- 
malization of the dark matter signal strongly depends 
on assumptions about the mass distribution in halos and 
subhalos, and is degenerate with the assumed value of the 
annihilation cross section. We have normalized the dark 
matter spectra in our scenarios to illustrate the potential 
of the anisotropy energy spectrum to identify a Galactic 
dark matter component while requiring that our models 
are both consistent with EGRET's measurement of the 
EGRB and plausible under standard assumptions 
about subhalo properties. 

For the extragalactic component of the diffuse isotropic 
emission, we consider unresolved blazars to be the only 
important contribution, and emphasize that this is the 
most relevant guaranteed source class since it is expected 
to dominate both the intensity of the EGRB above E > 1 
GeV and the angular power spectrum of the EGRB at 
I > 100. We obtain the intensity energy spectrum of un- 
resolved blazars as in [l9[ • Predictions for the amplitude 
of the blazar contribution to the isotropic diffuse back- 
ground vary by more than an order of magnitude (e.g. 
[201 ] and references therein), so we choose the normal- 
ization for our examples from within the range allowed 
by well-motivated predictions. The shape of the inten- 
sity energy spectrum Ie from a population of unresolved 
blazars with power-law spectra can be parameterized by 
the mean spectral slope «o an d the width of the spectral 
index distribution cto . We choose these parameters to be 
consistent with constraints on the spectral index distri- 
bution derived by [2lj based on EGRET observations. 
Our "reference blazar model" has parameters ctQ = 2.35 
and (To = 0.15. To approximate the spectral feature due 
to EBL absorption, we use the analytic expression given 
in [22| which applies an energy- and redshift-dependent 
cut-off to the extragalactic source spectrum, and assume 
for simplicity that all blazars reside at a redshift zq. 

The anisotropy energy spectrum. — Wc define the 
anisotropy energy spectrum as the value of the angular 
power spectrum at a fixed multipole I as a function of 
energy; we fix £ = 100 in the example scenarios pre- 
sented. We consider the angular power spectrum Ce of 
intensity fluctuations SI(ip) = — (I)) /(I), where 

I(ip) is the intensity in the direction ip. The angular 
power spectrum is given by Ci = { \ai m \ 2 ), where ag m are 
determined by expanding 5I(?p) hi spherical harmonics, 
SI(tp) = J2e, m a imYi m (ip) . 

For Galactic dark matter substructure we use the an- 
gular power spectrum for an anti-biased subhalo distribu- 
tion with mass function slope a = 0.9 and M min = 10 M Q 
calculated using the procedure of Q. We estimate the 
value of the angular power spectrum of emission from un- 



resolved blazars at £ = 100 from Fig. 4 of [12(. In general, 
the predicted amplitude of the angular power spectrum 
from substructure is much greater than that from blazars. 
This can be understood by noting that both source pop- 
ulations consist primarily of point sources, which results 
in a noise-like angular power spectrum, the amplitude of 
which is inversely proportional to the number density of 
sources per solid angle. As a cosmological source class, 
blazars have a much higher number density than Galactic 
substructure, and thus produce less angular power. 

The total angular power spectrum from extragalactic 
sources Cf G and Galactic dark matter substructure CP M 
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where /eg and /dm are the (energy-dependent) frac- 
tions of the total emission from extragalactic sources and 
Galactic dark matter substructure, respectively. Since 
extragalactic sources and Galactic dark matter sub- 
structure are uncorrelated, the cross-correlation term 
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measured angular power spectrum is given by 
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where Cf is the angular power spectrum of the signal 
(here, C\ at ) and Wi = exp(—£ 2 a£/2) is the window func- 
tion of a Gaussian beam of width ov We take the frac- 
tion of the sky observed / s k y = 0.75, and multipole bins 
of A£ — 100. The noise power spectrum Cn is the sum 
of the Poisson noise of the signal and the background, 
C" N = (47r/ sky /iV s )(l + (JV b /JV 8 )), where N s and N h are 
the number of signal and background photons observed. 
In this analysis, the EGRB and gamma rays from Galac- 
tic dark matter substructure together constitute the sig- 
nal, while any other observed gamma rays, such as those 
from Galactic sources other than dark matter substruc- 
ture, are considered background photons. To assess the 
impact of plausible backgrounds on the measurement er- 
rors, we assume N-b/N s = 10. This ratio is a rough but 
generous estimate of the factor by which the Galactic dif- 
fuse emission predicted by the GALPROP conventional 
model exceeds our isotropic diffuse intensity for latitudes 
|6| > 20° for energies up to - 50 GeV 0. We take the 
field of view of Fermi to be 2.4 sr, and assume an all- 
sky observation time of 5 years. The effective area A e f{ 
and angular resolution db of Fermi vary over the energy 
range considered here, so we approximate the energy de- 
pendence of these par ameters from the estimated Fermi 
performance plots 27 1. 

Fig. [l] shows an example intensity energy spectrum and 
corresponding anisotropy energy spectrum for m x = 700 
GeV, a mass motivated by attempts to explain the re- 
cent results of PAMELA and ATIC with dark matter 
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annihilation 2J, |25[ . The observed intensity energy spec- 
trum is the sum of the EBL-attenuated reference blazar 
spectrum and the dark matter spectrum, and in this ex- 
ample dark matter dominates the intensity energy spec- 
trum above ~ 20 GeV. The observed intensity spectrum 
is, however, also consistent with a blazar-only spectrum 
with a broader spectral index distribution (an "alterna- 
tive blazar model", ao — 2.32, ctq = 0.26) that has suf- 
fered EBL attenuation. In light of uncertainties in the 
properties of blazars, the EBL, and dark matter, the in- 
tensity energy spectrum alone is not sufficient to distin- 
guish between these two possibilities. In this case the 
anisotropy energy spectrum can break the degeneracy: if 
unresolved blazars were the sole source of the isotropic 
diffuse emission, the anisotropy energy spectrum would 
be constant in energy, but the presence of a dark matter 
contribution that varies with energy results in a modula- 
tion of the anisotropy energy spectrum. 

Fig. [2] presents a scenario with m x = 80 GeV, which is 
generally considered a more favorable mass for detection 
by Fermi. However, in this scenario the dark matter in- 
tensity is always subdominant, and as before the observed 
cut-off in the intensity energy spectrum occurs at an en- 
ergy consistent with EBL suppression of the EGRB, pro- 
ducing an acute degeneracy between the reference blazar 
model plus a dark matter contribution and an alterna- 
tive blazar model (ao = 2.28, Co = 0.26) without dark 
matter. Again, the anisotropy energy spectrum provides 
a means of robustly identifying a dark matter contri- 
bution: even though Galactic dark matter substructure 
never dominates the intensity energy spectrum, it pro- 
duces a strong feature in the anisotropy energy spectrum. 

In both examples, the error bars become prohibitively 
large for E < 1 GeV due to the angular resolution of 
Fermi below this energy, and at sufficiently high ener- 
gies due to lack of photons. In between these two regimes, 
the noise term in Eq. [2] (C-^/Wf) is negligible, and the 
uncertainties are quite small, governed primarily by the 
sample variance at the selected multipole. As a result, 
the departure of the measured anisotropy energy spec- 
trum from an energy-invariant quantity can be identified 
with high confidence, clearly indicating a transition in 
energy between source populations. 

We comment that the blazar intensity spectra (as well 
as the total intensity) in our examples fall noticeably be- 
low the EGRET data points. This reflects the expecta- 
tion that Fermi, with its enhanced point-source sensitiv- 
ity, will resolve a large number of extragalactic sources 
that had contributed to EGRET 's measurement of the 
EGRB, and consequently will measure a lower amplitude 
diffuse background. The EGRET data points are plotted 
to explicitly demonstrate that our models do not violate 
existing constraints. 

Discussion. — The observation of a modulation in the 
anisotropy energy spectrum robustly indicates a change 
with energy in the spatial distribution of contributing 




FIG. 1: Top panel: Example measured isotropic diffuse inten- 
sity spectrum. Shown individually are the spectra of Galactic 
dark matter substructure for m x = 700 GeV, the reference 
blazar model without and with EBL attenuation (zq = 0.4), 
and the unattenuated alternative blazar model. The 'total' 
signal is the sum of the attenuated reference blazar spectrum 
and the dark matter spectrum. Fhe EGRET measurement 
of the EGRB is plotted for reference (black crosses). Bottom 
panel: Energy dependence of the angular power spectrum of 
the total isotropic emission at multipole i = 100 for the sce- 
nario shown in the top panel. The anisotropy energy spectrum 
of Galactic dark matter substructure, unresolved blazars, and 
the total signal are shown. 




FIG. 2: Same as Fig.[T] for m x — 80 GeV. For the attenuated 
reference blazar model, zq = 1. 



source population(s). Although we have considered only 
the contributions of Galactic dark matter substructure 
and unresolved blazars to the isotropic diffuse back- 
ground, sources other than those explicitly considered 
here (e.g., Fermi irreducible backgrounds, the smooth 
dark matter halo, and additional extragalactic popu- 
lations including dark matter) which could induce an 
energy-dependence in the total angular power spectrum 
are not expected to provide significant power at the angu- 
lar scales of interest. Here we have not explicitly consid- 



4 



ered the possibility of energy dependence of the angular 
power spectrum of a single source class, although such 
a dependence could be produced by redshifting of hard 
features (e.g., lines, exponential cut-offs) in cosmologi- 
cal source spectra [ll|, [26| , spatial dependence or redshift 
evolution of the properties of source populations, or EBL 
attenuation. These mechanisms affect the anisotropy en- 
ergy spectrum by changing the relative contribution of 
different members of a given source class as a function of 
energy. For example, EBL attenuation of high-redshift 
blazars could reduce the number density of contribut- 
ing sources at sufficiently high energies. However, all of 
these mechanisms for generating strong features in the 
anisotropy energy spectrum would be accompanied by 
features (e.g., suppression, in the case of EBL attenua- 
tion) in the intensity energy spectrum at the same en- 
ergies, so these scenarios can be verified or excluded by 
joint examination of the intensity and anisotropy energy 
spectra. In general, confident identification of a specific 
source population will require information from the in- 
tensity energy spectrum since the modulation must have 
an energy dependence consistent with expectations for 
that source population and with features in the inten- 
sity energy spectrum. Moreover, although the noise-like 
shape of the angular power spectrum from substructure 
is not unique to this source class, the combination of 
its shape and amplitude, in addition to the energy spec- 
trum and large-scale angular distribution of substructure 
(which lacks a strong correlation with Galactic struc- 
tures), is not reproduced by any other known or predicted 
source class. 

A non-negligible contribution to the isotropic diffuse 
emission from Galactic dark matter substructure is guar- 
anteed to produce a modulation of the anisotropy energy 
spectrum. The precise shape, amplitude, and prominence 
of such a feature are governed by the angular power spec- 
tra of the EGRB and Galactic dark matter and their 
contributions to the total intensity, which are subject to 
considerable uncertainty and thus may differ significantly 
from the values adopted here, and therefore the scenar- 
ios presented should be regarded as examples, not pre- 
dictions. However, we emphasize that any modulation in 
the anisotropy energy spectrum can be used to extract in- 
formation about contributing source populations: a firm 
prediction for such a feature is not required to infer the 
presence of multiple or changing populations. 

It has not escaped our notice that in the limit where 
the baseline extragalactic anisotropy level is observable, 
Eq. Q] provides a novel way of extracting the shape of the 
dark matter annihilation spectrum even in cases where 
the dark matter signal cannot be disentangled from other 
contributions in the observed intensity spectrum. 
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